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D , Abstract 

O . A measurement of the first four moments of the hadronic mass distribution in i? — > Xcii^ decays 

is presented for minimum lepton momenta varying between 0.9 and 1.6 GeV, using data recorded 
with the BABAR detector. Furthermore, a measurement of the inclusive electron energy spectrum 
^ . for semileptonic B decays together with a measurement of its first, second and third moments for 

H I minimum electron energies between 0.6 and 1.5 GeV is reported. We determine the inclusive B — > 

Xcii^ branching fraction, Bdu, the CKM matrix element \Vcb\, and other heavy-quark parameters 
from a simultaneous fit to the measured moments. 



Contributed to the Proceedings of the XXXIX*^ Rencontres de Moriond on QCD and High 

Energy Hadronic Interactions 
3/28/2004—4/4/2004, La Thuile, Italy 



Stanford Linear Accelerator Center, Stanford University, Stanford, CA 94309 
Work supported in part by Department of Energy contract DE-AC03-76SF00515. 



1 Introduction 

Moments of inclusive distributions and rates for semileptonic and rare B decays can be related via 
Operator Product Expansions (OPE) [1] to fundamental parameters of the Standard Model, such as 
the Cabibbo-Kobayashi-Maskawa matrix elements \Vcb\ and \Vub\ [2] and the heavy quark masses rrib 
and rric- These expansions in l/mh and the strong coupling constant ag involve non-perturbative 
quantities that can be extracted from moments of inclusive distributions. In the kinetic-mass 
scheme [3] for example, these expansions to order 0{l/ml) contain six parameters: the running 
kinetic masses of the b— and c— quarks, mb{iJ.) and mdlJ-), and four non-perturbative parameters. 
We determine these parameters from a fit to the moments of the hadronic mass and electron energy 
distributions in semileptonic B decays to charm particles, B -^ Xciv. These measurements are 
based on data recorded with the BABAR detector [4] at the T{AS) resonance. 

2 Measurement of Hadronic Mass Moments 

2.1 Event Selection and Simulation 

The measurement of hadronic mass moments is based on a sample of 89 million BB pairs corre- 
sponding to an integrated luminosity of 82 fb~^. The analysis uses T(4S') -^ BB events in which 
one of the B mesons decays to hadrons and is fully reconstructed (Breco) [5] and the semileptonic 
decay of the recoiling B meson (i?rccoii) is identified by the presence of an electron or muon. This 
approach allows for the determination of the momentum, charge, and flavor of the B mesons. 
Semileptonic decays are modeled by a parameterization of form factors for B -^ D*£~T' [6], and 
models for B -^ Di-u,D**£-V [7] and B -^ D'Kl'V,D*Tii.-V [8]. Monte Carlo (MC) simulations 
of the BABAR detector are based on GEANT4 [9] . 

Semileptonic decays are identified by the presence of exactly one electron or muon above a 
minimum cut-off energy -Ecut) measured in the rest frame of the -Brccoil meson recoiling against 
the -Breco- The hadronic system X in the decay B -^ Xlv is reconstructed from charged tracks 
and energy depositions in the calorimeter that are not associated with the -Breco candidate or the 
charged lepton. The neutrino four-momentum py is estimated from the missing four-momentum 
Pmiss = Pr(45) ~PBrcco ~Px—Pei where all momenta are measured in the laboratory frame. We impose 
the following criteria to ensure well reconstructed events: E^i^s > 0.5 GeV, |pmiss| > 0.5 GeV, and 
l-f'miss — IPmissll < 0.5 GcV. We sclcct 7114 signal events over a combinatorial background of 2102 
events. The mass of the hadronic system Mx is determined by a kinematic fit that imposes four- 
momentum conservation, the equality of the masses of the two B mesons, and constrains p'^ = 0. 

2.2 Extraction of Hadronic Mass Moments 

In order to extract unbiased moments (M^), we need to correct for effects that can distort the 
mass distributions. We use observed linear relationships between the measured (M^) and generated 
{M'^ *''"'^) values from MC simulations in bins of M^ ^™^ (see Fig. 1) to calibrate the measurement 
of M'^ on an event-by-event basis. Since any radiative photon is included in the measured hadron 
mass and our definition of Mx does not include these photons, we employ PHOTOS [10] to simulate 
QED radiative effects and correct for their impact (less than 5%) on the moments as part of the 
calibration procedure. To verify this procedure, we apply the calibration to the measured masses 
for individual hadronic states in simulated B -^ XcH-v decays, and compare their calibrated mass 
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Figure 1: Results of the {Mx) calibration procedure. The calibration data and fit results are shown 
by the lower dashed line (circles), the verification by the upper solid line (triangles). 

moments to the true mass moments. The result of this test is also shown in Fig. 1 for Mx, indicating 
that the calibration reproduces the true moments over the full mass range. Corresponding curves 
are obtained for M^, M^, and M^. We observe no significant mass bias after calibration. The MC- 
based calibration procedure has also been validated on a data sample of partially reconstructed 
D*^ — > D^TT^ decays [5] where the low-momentum tt"*" serves as a tag and allows to select an 
inclusive event sample for which the true hadronic mass is known. 

2.3 Results and Systematic Errors 

The hadronic mass moments {M"^) after background subtraction, calibration and correction for 
detection and selection efficiencies are shown in the upper half of Fig. 3. The full numerical 
results can be found in [11]. The four moments increase as -Ecut decreases due to the presence of 
higher mass charm states. The moment measurements are highly correlated. The statistical and 
systematic errors are of comparable size. The dominant systematic error sources stem from the 
precision of the modeling of the detector efficiency and particle reconstruction, the subtraction of 
the combinatorial background of the i?reco sample and remaining B background, and uncertainties 
in the modeling of the hadronic states. 



3 Measurement of Electron Energy Spectrum and its Moments 

3.1 Event Selection 

The measurement of the electron energy spectrum and its moments is based on a data sample of 
47.4 fb~^ recorded at the T{AS) resonance and 9.1 fb~^ at an energy 40 MeV below the resonance, 
measured in the electron-positron center of mass frame. The analysis is similar to the BABAR 
measurement of the semileptonic branching fraction [12], but supersedes it by an order of magnitude 
in integrated luminosity. 

We identify BB events by observing an electron, etag, with a momentum of 1.4 < p* < 2.3 GeV/c 
in the T{AS) rest frame. These electrons make up the tag-sample which is used as normalization 
for the branching fraction. A second electron, e^ig, for which we require p* > 0.5GeV/c is assigned 
either to the unlike-sign sample (charge Q{etag) = —Q{e-sig)) or to the like-sign sample {Q{etag) = 
Qi^sig))- In events without B^B^ mixing, primary electrons from semileptonic B decays belong to 
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Figure 2: Electron momentum spectrum from B — > Xei'{'^) decays in the T{AS) frame after 
correction for efficiencies and bremsstrahlung, with combined statistical and systematic errors. 



the unlike-sign sample while secondary electrons contribute to the like-sign sample. We select BB 
events by making requirements on the charged and neutral multiplicities and event shape variables. 
Electrons originating from the same B meson as the tag electron are rejected by requiring cos a* > 
1.0 — p*(GeV/c) andcosa* > —0.2, where a* is the opening angle between the two electrons. 
Backgrounds from J/tp — > e~^e~ decays to the tag-sample, are suppressed by requiring the invariant 
mass Mee of the tag electron, paired with any electron of opposite charge with cos a* < —0.2, to 
be outside the interval 2.9 < Mge < 3.15GeV/c^. The efficiencies of these selection criteria are 
estimated by MC simulation. 



3.2 Measurement of the Electron Momentum Spectrum 

In order to extract the electron momentum spectrum remaining backgrounds have to be subtracted. 
Continuum background is subtracted from the tagged, like- and unlike-sign samples by scaling the 
off-resonance yields by the ratio of on- to off-resonance integrated luminosities, corrected for the 
energy dependence of the continuum cross section. Background electron spectra from photon 
conversions and Dalitz decays are extracted from data. Further backgrounds arise from decays of r 
leptons, charmed mesons produced in 6 ^ ccs decays and J/-0 or ip{2S) — > e~^e~ decays with only 
one detected e. These backgrounds are irreducible, and their contributions to the three electron 
samples are estimated from MC simulations. 

To account for B^B^ mixing, we determine the number of primary electrons in the i-th. p* bin 
from the like-sign and unlike-sign pairs as 
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where xo is the B'^B" mixing parameter xo = 0.186 ± 0.004 [13] and /o = B(T{4S) -^ B^B^ ^ 
0.490 =b 0.018 [13]. The parameter ej^ is the efficiency of the additional requirement on the opening 
angle for the unlike-sign sample. The spectrum of primary electrons is shown in Fig. 2 and is 
corrected for bremsstrahlung in the detector using MC simulation. 
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Figure 3: The measured hadronic-mass (a-d) and electron-energy (e-h) moments as a function of 
the cut-off energy, Ecut , compared with the result of the simultaneous fit (line) , with the theoretical 
uncertainties indicated as shaded bands. The solid data points mark the measurements included 
in the fit. The vertical bars indicate the experimental errors; i.e., the statistical and systematic 
errors added in quadrature; in some cases they are comparable in size to the data points. Moment 
measurements for different Ecut are highly correlated. 

3.3 Extraction of Electron Energy Moments 

Before extraction of the moments as a function of the minimal electron energy ranging from 
0.6 GeV to 1.5 GeV, the measured electron momentum spectrum has to be corrected for the 
contribution from B — > XuGV' decays. Defining Ri(Ecnt, fJ-) as J^ {Eg — fiy{dT/dEc) dEc, we 
present measurements of the partial branching fraction Mo(-E'cut)j the first moment Mi{Ecut) = 
Ri{Ecnt,0)/Ro{Ecnt,0) and the central moments M„(£'cut) = RniEcut,Mi{Ecut))/Ro{Ecnt,0) for n 
= 2, 3. Furthermore, the moments are corrected for the movement of the B mesons in the center of 
mass frame and QED radiative effects. The results are shown in the lower part of Fig. 3. The full 
numerical results can be found in [14]. The main systematic errors stem from uncertainties in the 
branching fractions of the irreducible background, electron identification efficiency and subtraction 
of the B -^ Xueu background. 



4 Determination of the Branching Fraction for B -^ Xclu Decays 
and of \Vcb\ from Hadronic Mass and Lepton Energy Moments 

The E'cut-dependent moment measurements discussed in the previous two sections are used to ex- 
tract the total branching fraction Bdv, |V^b| and other heavy quark parameters from a simultaneous 
X^ fit of OPE calculations in the kinetic mass scheme. In this scheme, the rate oi B ^ Xclv can 



Table 1: Fit results and error contributions from the moment measurements, approximations to 
the HQEs, and additional theoretical uncertainti es from Ug terms and other perturbative and 
non-perturbative terms contributing to Tdy 





\Vcb\{W-') 


m6(GeV) 


niciGcV) 


tiiiGeV') 


pUGeV') 


^G(GeV^) 


pisiGeV') 


BciA%) 


Results 


41.390 


4.611 


1.175 


0.447 


0.195 


0.267 


-0.085 


10.611 


Sexp 

Shqe 
5r 


0.437 
0.398 
0.150 
0.620 


0.052 
0.041 
0.015 


0.072 
0.056 
0.015 


0.035 
0.038 
0.010 


0.023 
0.018 
0.004 


0.055 
0.033 
0.018 


0.038 
0.072 
0.010 


0.163 
0.063 
0.000 


Stot 


0.870 


0.068 


0.092 


0.053 


0.029 


0.067 


0.082 


0.175 



be expressed as [15] 
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to 0{l/m^). The leading non-perturbative effects arise at 0{l/m'^) and are parameterized by ^J^^^i) 
and fiQ^fi), the expectation values of the kinetic and chromomagnetic dimension-five operators. 
At 0{l/ml), two additional parameters enter, p\){jx) and p\s{l^)-, the expectation values of the 
Darwin {D) and spin-orbit {LS) dimension-six operators. These parameters depend on the scale 
/i that separates short-distance from long-distance QCD effects; the calculations are performed 
for /i = 1 GeV [3]. The ratio r = m^/m^ enters in the phase-space factor ZQ{r) and the function 
d{r) ^ . HQEs in terms of the same heavy-quark parameters are available for the hadronic mass and 
electron energy moments [17]. Since many of these individual moments are highly correlated we 
select for the fitting procedure a set of moments for which the correlation coefficients are less than 
95%. Thus we only use half of the 28 mass moments, and retain 13 of the 20 energy moments. 

The global fit takes into account the statistical and systematic errors and correlations of the 
individual measurements, as well as the uncertainties of the expressions for the individual moments. 
The resulting fit, shown in Fig. 3, describes the data well with x^ = 15.0 for 20 degrees of freedom. 
Table 1 lists the fitted parameters and their errors ^. An additional error on [V^^l has been derived 
from the limited knowledge of the OPE expression for the decay rate, including various perturbative 
corrections and higher-order non-perturbative corrections [15]. As expected, mf, and rric are highly 
correlated and for the mass difference we obtain m\, — rric = (3.436 ± 0.025exp i Q.^IShqe i 
O.OlOoJGeV. 

Several crosschecks have been carried out to ensure that the fit results are unbiased. These 
include variations of the theoretical uncertainties and the set of moment measurements used in the 
fit. In particular these have been split into sets above and below a cut-off energy of 1.2 GeV and 
hadron mass and lepton energy moments. All results agree with each other within errors. Figure 4 
shows the Ax^ = 1 ellipses for |V^f,| versus rrn, and ;U^, for a fit to all moments and separate fits to 
the electron energy moments and the hadronic mass moments, but including the partial branching 
fractions in both. 



^Analytical expressions for zo{r) and d{r) can be found in [16]. 
^The correlations between the fit parameters can be found in [16]. 
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Figure 4: Fit results (crosses) with contours corresponding to A^^ = 1 for two pairs of the eight 
free parameters a) nih and b) //^ versus |V"cfe|, separately for fits using the hadronic -mass, the 
electron-energy, and all moments. 



5 Conclusions 

From the measured hadronic mass and electron energy moments we have determined the semilep- 
tonic branching fraction, \Vcb\ and the heavy quark masses rrih and ttt-c- 

iKfel = (41.4 ± 0.4e,p ± 0.%Qs ± 0.6iO X 10-3, 

Bceu = (10.61 ±0.16e:.p±0.06HQs)%, 

m5(lGeV) = (4.61 ± 0.05exp ± 0.04HQi? ± 0.02*,,) GeV, 
m,(lGeV) = (1.18 ±0.07e:..p±0.06HQi?±0.02t/,) GeV, 

In addition, the non-perturbative parameters in the kinetic scheme were determined up to order 
1/m^ without applying any external constraints and the results are in agreement with theoretical 
estimates. Consistent results were found from separate fits using hadronic mass or lepton energy 
moments only which gives confidence in the reliability of the OPE calculations. 
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